Quantum refrigerator driven by current noise 
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We proposed a scheme to implement a self-contained quantum refrigerator system composed of 
three rf-SQUID qubits, or rather, flux-biased phase qubits. The three qubits play the roles of the 
target, the refrigerator and the heat engine respectively. We provide different effective temperatures 
for the three qubits, by imposing external current noises of different strengths. The differences of 
effective temperatures give rise to the flow of free energy and that drives the refrigerator system 
to cool down the target. We also show that the efficiency of the system approaches the Carnot 
efficiency. 
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Introduction. — It is an interesting problem to discuss 
how small we can create a cooling machine and what 
would happen when quantum effects are taken into con- 
sideration, e.g., whether a quantum refrigerator could ex- 
ceed the classical Carnot efficiency. In practice, it is also 
a great challenge how to obtain lower temperature to im- 
plement quantum tasks. 

A lot of work has been done, both theoretically and 
experimentally [J-Q, to cool a quantum system. Most 
of the proposals require external control of microwave 
field for excitation or periodic control, and the cooling 
efficiency in experiments is usually quite low. 

However, recently, Linden et al. proposed a self- 
contained refrigerator system, which has a heat engine 
inside to drive the whole system The system con- 
tains three qubits, which play the roles of the target to 
be cooled, the refrigerator and the heat engine, denoted 
by 1, 2 and 3 respectively. The three qubits work un- 
der different temperature conditions, and they interact 
through 

V = Vg(|010)(101| + 1 101) (010|). 

Ei + E3 = E2 is also required to guarantee the energy 
conservation, where E a is the energy of each two-level 
system, H a = £' Q |1) ? (1|. 

The first term in V describes the main cooling process 
of qubit 1 (the target) by qubit 2 (the refrigerator), while 
the second term describe the reversal effect. In order to 
make sure the cooling process dominates, we must guar- 
antee that the third qubit (the heat engine) is hot enough. 
That is to say, the population of |1) of qubit 3 must be 
large enough, so as to enhance the cooling process. They 
analyzed the problem in a phenomenological dissipation 
model and demonstrated the cooling effect. Later works 
show that the efficiency approaches the Carnot engine 

In this paper, we propose a realizable model to imple- 
ment this self-contained quantum refrigerator by three 
rf-SQUID qubits, or rather, flux-biased phase qubits. 
Josephson circuits techniques are relatively sophisticated 



nowadays, no matter the fabrication, control or measure- 
ment, especially when we are dealing with quantum tasks 
that includes not too many qubits [lMil- Moreover. 
Josephson circuits have some unique merits to build this 
self-contained refrigerator system. 

One of the obstacles to build this refrigerator is that, in 
experiments, it is hard to build systems with 3-body in- 
teraction directly. However, indirect 3-body interaction 
may arise from the transmission of basic 2-body interac- 
tions, if we impose proper detuning condition. This is 
a method frequently used in optical lattice systems fl3l | . 
Fortunately, we can achieve such indirect 3-body inter- 
action in our system proposed in this paper. 

Besides, another problem that seems more difficult is 
how to maintain the three microscopic qubits in different 
temperatures, which separate from each other usually at 
a distance of only several microns. We would solve it 
by utilizing effective temperatures of the Nyquist noises 
brought in from external circuits. 

In Josephson circuits, one of the important sources of 
noise comes from external currents. The current noises 
are unavoidably brought in from external circuits that 
perform controlling or measuring, and transfer along the 
circuit wires. Usually, we suppress current noise by filters 
or delicate design to enhance the quantum coherence of 
Josephson qubits. 

In our proposal here, we pour into the three qubits with 
noises of different strengths. Current noises run along 
wires from outside, exactly as hot water runs along pipes 
from a heat source. Therefore, equivalently, we offer the 
three qubits with independent different thermal reser- 
voirs. Moreover, the equivalent thermal reservoirs that 
we provide here are stable, whose temperatures cannot 
be affected by the three qubits. In experiments, Joseph- 
son circuits are usually thermally anchored to the mix- 
ing chamber in a dilution refrigerator at a temperature 
T m i x ~ 10 mK, while the effective temperature caused 
by the current noise may be even as high as ~ 300 mK 
[14j . We believe that we can make use of such effective 
temperature of noises in our self-contained refrigerator 
system. 
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Figure 1: The design of circuit of our flux- biased phase qubits 
system. The part inside the dashed line square represent three 
identical mutual inductance coils interacting with each other 
symmetrically. External inductances provide flux biases for 
each qubit. 



Proposal. — We show our circuit design in Fig. 1. Flux 
bias for each qubit is provided by external noisy current, 
and we assume this is the main contribution to dissi- 
pation in our system. The three qubits interact through 
mutual inductive coils overlayed together (ill . [l2| , as rep- 
resented in the dashed line squre part in Fig. 1. For sim- 
plicity, we assume the three mutual inductive coils are 
identical. 

To get the quantum description of the system, we can 
write down the classical equations of motion of the cir- 
cuit, then get the Lagrangian of the system, and finally, 
we can obtain the conjugate canonical momentum and 
the Hamiltonian, 



a = l 



ZUq /i±Jry 



$2(L M + M) 



(Lm + 2M)(L m -M) 
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a=l a<f) 



I) 



where $o = 3>o/27r and $o is the flux quantum. $° xt is 
the external flux imposed to the rf-SQUID loop. Q a is 
the charge carried by the capacitance of the Josephson 
junction, and ip a is the superconducting phase difference 
across each junction. L a is the self-inductance of each rf- 
SQUID loop. Lm and M are the self and mutual induc- 
tances of the mutual inductive coils overlayed together, 
as denoted in the dashed lined square in Fig. 1. 

We set the parameters in such a way that each qubit 
works in a meta-stable cubic well approximately, and 
choose the lowest two bound state as a two-level system. 
Move the origin of ip a to the stable point of potential 
U a {(p a )i we can rewrite the Hamiltonian as 



H 



+ Vij, 

a=l 
3 

E 



m a uj n x n 



2m c 



g ^2 x &xp(2) 



m a , w a , A and the coupling constant g are determined 
from the parameters in Eq. (pj. x a = <p a — <^ ta " is the 
translated coordinate whose stable point is settled at the 
origin point . 

We can treat H a as a harmonic oscillator plus a 
cubic term as the perturbation. So we have x a = 
[h/2m a uj a }2 (a a + aJJ, and a a is the annihilation oper- 
ator of the harmonic oscillator. Modifications should be 
made to the energy and eigenstate of the corresponding 
harmonic oscillator. 

If we focus on the dynamics of the lowest two lev- 
els, we can write down the low-energy effective Hamil- 
tonian of each single qubit as H a = £' Q |1) Q (1|, where 
E a ~ 0.95huj a is the anharmonically modified energy of 
the qubit (More details of the computation can be found 
in the appendix in Ref. [11]) . 

Indirect 3-body interaction. — Although Vi n t in the 
Hamiltonian Eq. ^ only involves two-body interactions, 
we can still consider indirect interactions of higher order 
in interaction picture, and obtain the effective Hamilto- 
nian through 



Wint(i) = Texp 



ft Jo 



<1tHi(t) 



= l + -H cS t + o{t 2 ). 
n 

The elements of the effective Hamiltonian up to the sec- 
ond order is 



H. 



(n\V iD 



v int J 



-> tr( 2 ) 
n -n^ff m 



E k -E 



where E n = E m = E and \n) = \n1n2n3} is the state of 
the three qubits. 

We require that E\ + E3 — E2 and the three en- 
ergy are incommensurable. That implies that, besides 
the diagonal terms, the only transition process that re- 
mains is |010)(101| + h.c, when we apply rotating-wave- 
approximation (RWA), i.e., when we focus on the long- 
time effects. Since Vint involves only two-body inter- 
actions, we have (101 |Vj n t| 010) = 0. When dealing 
with the second order, we would come across terms like 
(n a \x a \n a }, which would not vanish when we add an- 
harmonic modifications to the states. Finally, we can get 
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an effective Hamiltonian, 

(3) 



a a</3 

+.g(|010)(101| + |101)(010| 



Here, we denote h a — |1) Q (1|. The coefficients of diago- 
nal terms is of order g, and the transition magni- 
tude g ~ g 2 . We can see that the interaction strength is 
quite small. 

Dissipation. — Dissipation due to external current 
noise has been well discussed in literature 15 -HJ. In 
our system, the flux bias is provided through <3?® xt = 
L® xt J a = L° xt [/ a + i a (t)]. L° xt comes from an external 
coil that provides flux bias for a single qubit, and i a (t) 
is the noise that satisfy (i a (t)) = 0. Tracing back to the 
single qubit term in the Hamiltonian of the Josephson 
circuit Eq. {TJ), we can figure out that the qubits inter- 
act with the current noise through an interaction term 
la<&a{t)ia{t)i where ^ a collects the parameters. 

The total Hamiltonian of the three qubits, each reser- 
voir and their interactions is 



Usb = n s + Hi 



rLh 



where H s = H eS , and H ia t = J2 a ^x a (t) ■ i a (t). Ub 
is usually treated as a collection of harmonic oscillators 
phenomenologically. 

We can write down the master equation of the sys- 
tem, after a complicated procedure of approximations 
that mainly include Born-Markovian approximation and 
RWA O G3|. The master equation of the thrce-qubit 
system is, 

d t p = -i[H cS ,p] + y^V a p 

(X 

= -i[H eS ,p] 

+ J2 T » ( N »{E a ) + 1) [2a a pal - {a a a a ,p}+\ 

a 

+ T a N a (E a ) [2a a pa a - {a a al, p}+] . (4) 

a 

Here, N a (E a ) = {exp((3 a E a ) — and a a are 

the correponding harmonic operators of each qubit. 
F a N a (E a ) oc Sf(E a /h), where Sf(uj) is the power spec- 
tral of current noise reflecting the dissipative impedance 
of the circuit, 

/>oo 

Sf(uj) = / e-**(i a (t)i a (p))dt. 
Jo 

Steady solution. — We concentrate on the equilibrium 
behaviour of the system here, and especially, we want to 
obtain the final steady distribution (h a ) of each single 
qubit. Then we can get the effective temperatures. 



We can get eight independent linear equations about 
{h a ), (hahp), (hih 2 h 3 ) and (Av) = g /aia|a 3 ^-h.c, by 
multiplying respective observables to the master equation 
and then tracing out the average. And we also have the 
steady solution, 

[Ua) ~ 2N a + l + 2r Q (2iV Q + l)' W 
where N a = N a (E a ). Let M a = T a (2N a + 1), we have 



(Av) = 



and 



ig 2 G 



X l +g 2 (X 2 +X 3 ) 
if (JViJV 3 -N 2 - NiN 2 - N 2 N 3 ), 



(NiN 3 -N 2 - NiN 2 - N 2 N 3 ) 



G = 4T 1 r 2 r 3 (M 1 +M 2 + M 3 )Y[(M a + M p ), 

a</3 

X 1 = 2(1 + ^-)M 1 M 2 M 3 (M 1 + M 2 + M 3 ) 2 
x H(M a + Mp), 

a</3 



Xo 



Xo 



AMiM 2 M 3 + M a Mp(M a + Mp) 

a<j3 

x *[[{M a + Mp), 

a<f3 

-T^M^iMi + M 2 ){M! +M 2 + 2M 3 ) 
+T 1 T 3 MiM 3 {Mi + Ma)(Mi + 2M 2 + M 3 ) 
-T 2 T 3 M 2 M 3 {M 2 + M 3 ){2Mi +M 2 + M 3 ), 

where B = Y ia M a , and A = D^> - + + D^ 13 \ 
The steady solution of (h a ) contains two terms. When 
there is no interaction between the three qubits, they 
would respectively decay into each Boltzmann distribu- 
tion, resulting from the weak coupling with their reser- 
voirs, as represented by the first term in Eq. ([5]). The 
second term is the contribution of interactions between 
qubits and that alter the populations. The effective tem- 
peratures would be changed as long as (Av) =^= 0. 

Cooling condition and efficiency. — When we say that 
we want to make the system run as a refrigerator in or- 
der to cool down qubit 1, we have actually implied that 
the initial temperature of qubit 1 is the lowest, other- 
wise, it could be cooled through heat transport. We can 
check that each term in £ is positive. Therefore, to lower 
the population number (hi), we must make sure that 
i (Av) < 0, and that gives the cooling condition, 



(ft ~ 02)Ek < (ft 



Ei ft 

E 3 ft - ft 



(3 3 )E 3 , 
ft 
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Remember that E\ + £3 = E%. As we know that T\ 
should be the lowest, so we must have T\ < T2 < T3. 

In order to compute the cooling efficiency of the refrig- 
erator, we have to compare the amounts of heat exchange 
of the target (qubit 1) and the heat engine (qubit 3) with 
their environments. We should go back to the master 
equation Eq. ((4]). The unitary term represents the con- 
tribution of doing works between the three qubits, and 
the second dissipative one represents the heat exchange 
with the environment. Therefore, we can get the heat 
exchange of each qubit with their environments per unit 
time by 

Q a = Tr [E a h a ■ V aP ] = {-l) a i (Av) E a . 

Now we arrive at an interesting result that the effi- 
ciency of the system is 

„q = 9e. = Ei < = „q 

1 Qh E 3 < - l - W 

This is exactly same with the results in previous work in 
Ref. 0. In their work, they also figure out that this is the 
upper bound on the efficiency of any such engine running 
between three reservoirs which extracts heat from the 
bath at Tq using a supply of heat from the bath at Tjj ■ 
We have to emphasize that T a here is the temperature 
of the reservoir where each qubit stays, and it is com- 
pletely determined by the strength of noise from external 
circuit. The change of the population of the qubits does 
not affect the temperature of reservoirs, which can be 
only controlled through external filter circuits. The ther- 
mal temperature of the qubits provided by the dilution 
refrigerator is the same, while the strengths of random 
motions of the electrons inside the conducting wires are 
different between the three qubits, and that give rise to 
an equivalent effect of providing different thermal condi- 
tions. 

In Ref. , it was mentioned that there is no other 
theoretical limitation for the refrigerator even when ap- 
proaching the absolute zero. However, in our system 
implemented by Josephson circuit proposed here, we 
have some additional constraints to the working condi- 
tions. First, the temperature of the heat engine Th must 
not be too high comparing with the excitation energy 
E3 ~ 10 GHz ~ 1 K, otherwise, population of higher en- 
ergy levels must be taken into account, and we cannot 
treat our system simply as a two-level system. Second, 
if we greatly suppress the external noise of qubit 1 in or- 
der to make Tc lower (~ 50 mK), the low frequency 1/f 
noise would dominate in the dissipation, therefore, our 
analysis above based on Markovian approximation in in- 
valid. After all, it has to be tested in realistic experiments 
how cold the refrigerator proposed here could achieve. 
However, this could still provide an astonishingly high 
efficiency when we compare it with macroscopic classical 
heat machines. 



Summary. — In this paper, we proposed a Josephson 
circuit system to implement a self-contained refrigerator. 
We demonstrated the configuration of the circuit, ana- 
lyzed the dynamics of the system and got the final steady 
distribution. By controlling the different strengths of 
the noises pouring into the qubits, surely the refrigera- 
tor would cool down the target with high efficiency even 
approaching up to the Carnot up-bound. We believe our 
proposal is realizable with present technology. 

The system contains a microscopic heat engine to 
maintain the running of the whole system automati- 
cally. The power actually comes from the external circuit 
noises. We do not need direct control of external control 
such as periodic microwave fields. What surprised us is 
that such a microscopic quantum heat machine driven by 
noises could run with quite high efficiency. This gives us 
an inspiration that maybe we can develop more quantum 
devices driven by noises. No need to say, there is a great 
advantage when there is a self-contained energy provider 
in microscopic quantum tasks. What's more, noises of 
many types may be utilized, such as noises from external 
current, flux or light, maybe even from inner fabrication 
defects. 

When our work was almost completed, we noted that 
Mari and Eisert proposed a similar mechanism of cooling 
by heating in an opto- mechanical system [l9j . 
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